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Scaling Relations of Star-Forming Regions: from kpc-size 
clumps to Hll regions 
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ABSTRACT 

We present the properties of 8 star-forming regions, or 'clumps,' in 3 galaxies at 
z ~ 1.3 from the WiggleZ Dark Energy Survey, which are resolved with the OSIRIS 
integral field spectrograph. Within turbulent discs, a ~ 90 km s" 1 , clumps are 
measured with average sizes of 1.5 kpc and average Jeans masses of 4.2 x 10 9 
M Q , in total accounting for 20 — 30% of the stellar mass of the discs. These findings 
lend observational support to models that predict larger clumps will form as a 
result of higher disc velocity dispersions driven-up by cosmological gas accretion. 
As a consequence of the changes in global environment, it may be predicted that 
star-forming regions at high redshift should not resemble star-forming regions locally. 
Yet despite the increased sizes and dispersions, clumps and Hll regions are found to 
follow tight scaling relations over the range z = — 2 for Ha size, velocity dispersion, 
luminosity, and mass when comparing > 2000 Hll regions locally and 30 clumps at 

T 4.18±0.21 and L K M 1.24±0.C^ 



z > 1 (a oc r u - 42±u ' U3 , L Ha ex r^ ±au4 , L Ha oc a*™*"^, and L Ha oc M^ UM0 ). 
We discuss these results in the context of the existing simulations of clump formation 
and evolution, with an emphasis on the processes that drive-up the turbulent motions 
in the interstellar medium. Our results indicate that while the turbulence of discs 
may have important implications for the size and luminosity of regions which form 
within them, the same processes govern their formation from high redshift to the 
current epoch. 

Key words: galaxies: formation - galaxies: high-redshift - galaxies: star clusters - 
galaxies: kinematics and dynamics - ISM: H II regions. 



1 INTRODUCTION 

Star-forming galaxies at high redshift are gene rally clumpy, 
gas-r i ch, and highly efficient at forming stars (jDaddi et all 
l20ld : iTacconi et all l201Ct ICombes et all l201lh . Compared 
to local star-forming galaxies that contain on average hun- 
dreds of Hll regions, galaxies in the early universe are made 
up of a handful o f kilo-parsec sized star-form ing 'clumps' 
jCowie et all 1 19951 : Ivan den Bergh et al.lll99ot ). These mor- 
phologies cannot be explai ned by simple re solution effects or 
an imaging K-correction (|Dickinsonl 120001 ) . Artificially red- 
shifting images of local star-forming galaxies with distribu- 
tions of hundreds of H II regions produces high-redshift im- 
ages with mostly even light distributions (|Elmegreen et all 
120091 ). Furthermore, clumpy structures have been seen in 
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broadband images as well as UV emission at high redshift, 
roviding evidence for the bu i ld-up of massive star clusters 
Forster Schreiber et alll2010l : lGuo et al.ll201ll ). 

New kinematic data from integral field spec- 
troscopy has supplemented the morphological studies 
to reveal that the clumps are embedded in massive 
rotating discs with high gas velocity dispersi ons of 



50 — 200 km s 1 (e.g . Ivan Starkenburg et all 
120091: iForster Schreiber et al 



Swinbank et al 



2008; 



2009: 



Jones et all l20ld : iGenzel et all l201ll : IWisnioski et all 
20111 ). The clumps are found to have comparable ve- 
locity dispersions to the global disc dispersion with 
typic al sizes of ~1 kpc and ma s ses of ~ 10 8 — 10 10 
M n (lElmegreen fc Elmegreenl 120051: lElmegreen et all 120091: 



Genzel et alll201ll : lForster Schreiber et allboilUGuo et all 
20111 ). 

The kinematic observations are in broad agreement 
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with models of unstable disc formation. In these models, 
motions in the turbulent interstellar medium (ISM) result 
in high local velocity dispersions that drive up the Jeans 
mass and Jeans length resulting in massive star-forming re- 
gions in gravitationally unstable discs. The timescale of the 
clumpy phase is predicted to be short, of order ~ 0.5 Gyr. 
The interactions of the clumps and high turbulence heats 
the gas resulting in a thick disc with scale height of ~ 1 
kpc. The more massive clumps then spiral to the centre 
of the galaxy due to dynamical fr ic tion to create a bulge 
llNoguchilll999l;llmmeli et al.ll2004aH bl; iBournaud et alj|2007t 
lElmegreen et al.ll2008l )~ 

In more recent models the turbulent discs are con- 
tinuously fed by cold gas accretion which allows multiple 
episodes of instability. In these simulations, the clumpy 
phase can have a lifetime of a few Gyr. Individual clumps 
last for ~ 0.5 Gyr if they have high enough star for- 
mation rates to survive destruction by stellar feedback 
(JKrumholz fc Dekell 120101 ) . Total star formation rates are 
approximately equal to the accretion rate with most of the 
galac tic star formation occurring in the clum ps (~ 100 M© 
vr" 1 ; iDekel et alj|2009l ; ICeverino et alj|2010h . The mass of 
individual kpc-size clumps, ~ 10 8 — 10 10 Mq, are a few per 
cent of the mass of t he disc amounting to a total of 20% of 
the mass of the disc (|Dekel et al.l 12009 ; iKrumhoIz fc Dekell 
l20ld : lGenel et alj|20ich ~~ 

However, there is disagreement amongst the models 
as to how clumps are able to survive long enough to be 
prevalent in observations. Clumps with supersonic veloc- 
ity dispersion and high star formation rates are expected 
to ex pend all their gas within a disc dynamical time (50 
Myr; IDekel et al.l 12009 ) unless something supports them 
from collapse. Turbulence, a natural candidate, is needed 
to maintain the high local velocity dispersions and thus the 
large instability scale. However, the restoring force of tur- 
bulence is uncertain with possible sources inc luding grav- 
itational motions of clumps w it hin the disc ( Dekel et al.l 
12009" ; iKrumholzfc Dekell _ 
self-gravity (|lmmeli et al.l 



201$., 
it all | 



disc 
Bournaud et all 120071 ; 



201 01; lAumer et al 
\2004\JT , 
ICeverino et all l20ld ; iBurkert et all |201Ct). radiative feed 
back JMurravet all l20ld; IKrumhoIz fc Dekell l20ld). su 



perno vae ( Mac Low fc Klessen) 2004 ; Ijoung fc Mac Low! 
2006) , and gas accretion ([Fqrster Schreiber et al.l 120061 ; 



Genzel et al.ll2006l ; lGenel et al.ll2010l ). However, turbulence 



is not the only mechanism for stabilising the clumps. 
If the regions become rotationally bound they may be 
supp orted against total free-fall collapse by their rota- 
tion (|Ceverino et alj|201ll ; lAgertz et alll200d ; lAumer et all 
l20ld ). 

Furthermore, it is unknown how long individual clumps 
survive and how they reach their end. In the extreme envi- 
ronments expected in unstable gaseous discs, radiation pres- 
sure and momentum-driven galactic winds are expected and 
may be strong enough to disrupt the clumps before they are 
able to virialise an d fall to the centre to form a bulge v ia 
dynamical friction (|Murrav et al.ll20ld ; lGenel et alj|20ld ). 

Kinematic measurements of clump motions and turbu- 
lence are therefore essential for understanding how high- 
redshift galaxies form. Measurements of velocity dispersion 
of clumps have become more available with the advent of in- 
tegral field spectroscopy (IFS) and adaptive optics (AO) sys- 
tems |Swinbank et aklfeood ; I Jones et alj|20ld : iGenzel et all 



120111 ). These can be used to test clump formation models 
and the driving forces behind their survival by using gas 
velocity dispersion as a proxy for turbulence. 

Recent IFS observations have suggested that high- 
redshift clumps may form differently from local HII regions 
as a result of a more turbulent environment with higher star- 
formation efficiencies, and gas densities. Clumps in intrinsi- 
cally L* and sub-L, galaxies at z > 2 have been found to be 
two orders of magnitude more luminous than HII regio ns at 
a given radius (|Swinbank et al.ll2009 ; Ijones et al.ll20ld ) , and 
their clump gas velocity dispersions lie in the r ange 60 — 100 
kms -1 , 2 — 6x greater than local HII regions JGenzel et al.l 
120111 ). However, these properties are comparable to those of 
the largest HII regions locally, commonly called 'giant' H 
II regions such as 30Doradus and II Zw 40 which have the 
same offset in star formation from 'normal'-size HII regions 
as high-redshift clumps. This offset in luminosity indicates 
that both classes have higher star formation surface densi- 
ties than regular HII regions suggesting that the most apt 
comparison at the different epochs is between clumps and 
giant HII regions. 

In addition to HII regions and giant HII regions, Ha 
results of high-redshift clumps have also been directly com- 
pared to local giant molecular clouds (GMCs; lMurrav et al.l 
120101 ). Some simulations are based on the assumption that 
giant clumps at hi gh redshift represent sing l e star-forming 
mole cular clouds (|Krumholz fc Dekell l20ld ; iMurrav et al.l 
120101 ). However, only a few clumps have been observed at 
high redshift from molec ular gas by observing the regions 
in highly lensed galaxies (|Swinbank et al.ll20ld ). It remains 
unclear what are the best analogues of high-redshift clumps 
at low redshift and whether the clumps form under the same 
conditions as GMCs or HII regions. 

One strategy to compare star-forming regions at 
high and low redshifts is to look at scaling relations 
to see if they evolve with the cosmic time. An inves- 
tigation of this kind is justified by the literature as 
extensive research has been dedicated to the relation- 
ships found between three key properties: luminosity, size, 
and velocity dispersion for giant HII regions and GMCs 
(e.g. iTerlevich fc Meln"icklll98ll; [Gallagher fc Hunter! Il983| ; 



Larsonl Il98ll; lArsenault fc RovMl988i; iRozas et al J 12009; 



Bastian et all 120061 ; iMonreal-Ibero et al.l 120071 ). The rela- 



tionships between size, luminosity, and velocity disper- 
sion have been greatly debated with proportionalities rang- 
ing from a oc r i-^ -3.68 L K ,.1.92-3 L K ^2.6-6.6 

JFuentes-Masip et alJlbood i lGutierrez et alj|201ll ). However, 
these discussions often suffer from small-number statistics 
and biased measurements. We will show that, by bringing 
together multiple studies from the literature, measured in a 
self-consistent way, the spectrum of HII regions form rela- 
tively consistent relationships between size, luminosity, and 
velocity dispersion, which can be extended to include high- 
redshift clumps. The scaling relationships observed have im- 
plications for how star-forming regions form in the changing 
cosmic environment. 

In this paper we present properties of eight clumps from 
three z ~ 1.3 galaxies in the WiggleZ survey . These galax- 
ies were selected from I Wisnioski et al.l ((2011), herein Paper 
I, in which we presented Ha kinematics of 13 star-forming 
galaxies observed with IFS. In Section 2 we introduce the 
WiggleZ clumps and a comparison sample at low and high 



© 2011 RAS, MNRAS 000, [TV?? 



Scaling Relations of Star-Forming Regions 3 



Table 1. High-Redshift Clump Samples 



Paper 




z 


Parent 


Number of 


Number of 


Radius 










Sample 


Clumps 


Galaxies 


Method" 


This Work 




1.28- 


-1.46 


WigglcZ 


8 


3 


core 


Jones et al. 2010 




1.68- 


-2.65 


lensed 


8 


1 


isophote 


Genzel et al. 2011 




2.18- 


-2.26 


SINS 


5 


3 


core 


Forster-Schreiber et al. 


2011 


2.26 


SINS 


7 


1 


photometric 


Swinbank et al. 2009 




4.92 


lensed 


2 


1 


isophote 



A more detailed description of how regions were measured is detailed in the Appendix. 



Table 2. Low-Redshift HII Region Samples 



Paper 


Number of Galaxies 


Parent Galaxy 


N a 


Radius 


Dispersion 




in Parent Sample 


Types 




Method 


Method 


SINGS (Kcnnicutt et al. 2003) 


7NGC 


Spirals, Irregulars 


2091 


core 




Gallagher & Hunter 1983 


7 NGC; Mk 35; 
Haro 22; A1004+10 


Irregulars 


30 


core 


echelle 


Arsenault et al. 1988 


26 NGC; 5 IC; 5 HO 


Spirals, Irregulars, 
Dwarf irregulars 


57 


isophotal 


Fabry-Perot 


Bastian et al. 2006 


NGC 4038/39 (Antennae) 


LIRG 


4 


isophotal 


IFS 


Rozas et al. 2006 


10 NGC 


Isolated spirals 


13 


40% isophotal 6 


echelle 


Monreal-Ibero et al. 2007 


5 IRAS 


ULIRGs 


12 


half-light 


IFS 



a Number of HII regions used in this paper. 

6 A method that measures core light. The size of the 40% isophote is equal to the diameter of the circle containing the same area as 

40% of the peak region surface brightness. 



redshift. We address the most reproducible measures of size, 
luminosity, and velocity dispersion for comparison across 12 
Gyrs (z = — 5) in Section 3 and apply them to the Wig- 
gleZ clumps. In Section 4 we quantify disc instabilities and 
estimate the mass of star-forming regions. In Section 5 we 
investigate the scaling relations of size, luminosity, and ve- 
locity dispersion of star-forming regions and their possible 
theoretical drivers. In Section 6 we use the results of Sec- 
tions 3—5 to test predictions from models and simulations 
of clump formation. In Section 7 we summarise the prop- 
erties of the WiggleZ clumps presented here. A standard 
ACDM cosmology of fi m = 0.27, fi A = 0.73, h = 0.7 is 
adopted throughout this paper. In this cosmology, at red- 
shift z — 1.3, 1 arcsec corresponds to 8.6 kpc in physical 
coordinates. 



2 DATA 

In this Section we present 8 new clumps, at z > 1, 
in g alaxies selected from t he WiggleZ Dark Energy Sur- 
vey (|Drinkwater et alj|20ich and bring together comparison 
samples of star-forming regions at high and low redshift. For 
the comparison samples we only consider clump and HII re- 
gion data from Hq observations for consistency and due to 
the wealth of data available in the literature. We focus on 
studies with published Ha fluxes, sizes, and velocity disper- 
sions for individual regions. A summary of the comparison 
data is given in Table [l] (high redshift) and Table [2] (low 
redshift), corrections to the distances, sizes, and cosmology 
that we introduce to improve the comparison of different 
samples are detailed in the Appendix. 



2.1 WiggleZ Clumps atz~ 1.3 

The new clump data presented in this paper are derived from 
the IFS observations of Ha luminous galaxies at z ~ 1.3 first 
introduced in Paper I. Data was taken with OSIRIS (OH 
Supp ressing InfraRed Imaging Spectrograph; lLarkin et al.l 
120061) with the Keck II Laser guide star adaptive o ptics sys- 
tem (|Wizinowich et al.ll200rj ; Ivan Dam et alj|2006h . OSIRIS 
is a lenslet array spectrograph with a 2048 x 2048 Hawaii- 
2 detector and spectral resolution R ~ 3600. All galaxies 
were observed with the 0.05" pixel scale. The FWHM of the 
tip-tilt stars ranged from one and a half to two pixels, 0.062- 
0.1", with an average Strehl ratio of 30% estimated from the 
tip-tilt stars. 

Of the 13 galaxies presented in Paper I, three galaxies 
have Ha morphologies with multiple resolved star-forming 
regions. A summary of the global properties of these galaxies 
are given in Table [3] From this data we measure individual 
Ha clumps and their kinematic properties. Clumps are iden- 
tified solely from 2D Ha emission images. The unsmoothed 
raw Ha images used throughout this analysis were created 
from a summation over the spectral dimension within 1- 
sigma of the systemic velocity, as determined by a single 
Gaussian fit to the integrated spectrum. None of the galax- 
ies have detectable continuum that could contaminate the 
Ha images. We find 8 resolved clumps in total. The identi- 
fication of clumps and measurement of their properties are 
discussed in Section 3. 



2.2 High-redshift comparison sample 

At high redshift the study of clumps suffers from low-number 
statistics due to the resolution limits of current instrumenta- 
tion from the ground. The high-redshift comparison sample 
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Table 3. Global Properties of WiggleZ Clumpy Galaxies 



ID 


z 


fa 

J Ha 


log[M„(M )] 


^mean 
(km s- 1 ) 


""net 

(km s- 1 ) 


''•'shear 

(km s- 1 ) 


12 + log(0/H) 


N" 


WK0912.13R 
WK1002.61S 
WK0909.06R 


1.2873 
1.3039 
1.4602 


48 ± 16.1 
81.5 ± 5.3 
94.9 ± 6.7 


10.7 ± 0.2 
10.7 ± 0.7 
11.0 ± 0.2 


81.6 ± 27.7 
85.1 ± 20.3 
92.6 ± 27.8 


104.3 ± 27.1 
88.4 ± 4.4 
153.3 ± 8.4 


130.5 ± 11.0 
71.5 ± 10.7 
160.8 ± 2.4 


<8.47 

<8.37 

8.55 ± 0.1 


2 
3 
3 



"Emission line flux in units of 10 ergs s 1 cm 
''Number of clumps measured in each galaxy. 



consists of 22 clumps at z > 1.6 as given in Table [T] which 
are all at higher redshift than the WiggleZ clump s . The 
data for two regions are taken from ISwinbank et al.l (| 20091 ) 
for the lensed z — 4. 92 galaxy C 1 1358+ 62. Eight regions 
are obtained from the I Jones et al.l (|20T0t ) lensed sample at 
z ~ 2, and twelve r egions are taken from the SI NS survey 
JGenzel et alj|201ll ; lF5rster Schreiber et al.ll201ll ) at z ~ 2. 
We also compare to simulated high-redshift (z ~ 1.9 — 3.0) 
clumps extracted from five simulated gravitationally unsta- 
ble discs with baryonic masses of 10 10 — 2 x 10 11 M ^ at 
different stages of their lifetime (|Ceverino et al.ll201rh . 



3 REGION PARAMETERS 

To accurately compare the data from the WiggleZ clumps 
with the two comparison samples the clump properties from 
each sample should be measured in a consistent way. In this 
section we investigate the best methods to determine the 
radius, luminosity, velocity dispersion, and metallicity, and 
apply them to local HII regions in the SINGS galaxies and 
toz~ 1.3 clumps in the WiggleZ galaxies. We consider the 
effects of resolution on these measurements and the impli- 
cations they will have for the scaling relations studied in 
Section 5. Results are summarised in Table [4] 



2.3 Low- redshift comparison sample 

At low redshift we select a sample of extragalactic H 
II regions and a sample of giant extragalactic HII re- 
gions. The difference between 'normal' HII regions and 'gi- 
ant' HII regions is determined not only by greater size 
and luminosity but primarily by having supersonic turbu- 
lence (a > c BOU nd), a property that was later integrated 
into t he definition of giant HII regions (ISmith fc Weedmanl 
Il97d : iTerlevich fc Melnick|[l98H ; iGallagher fc Hunterlll983F 
The sound velocity for HII regions is defined as chii = 
y^kfimTun = 13 km s _1 where k is the Boltzmann con- 
stant, (j, is the mean molecular weight, m is the mass of hy- 
dr ogen, and Tun is the temp erature for HII regions (~ 10 4 
K: iFuentes-Masip et al.ll2000l ). For the velocity dispersions 
to be supersonic a mechanism is required to maintain the 
turbulence, which is expected to dissipate by shocks within 
a free-fall time, making them interesting objects to study. 

The comparison data for the normal HII regions are 
derived from the Spitzer Infra red Nearby Galaxies Survey 
(SINGS . iKennicutt et "aL1l2003h . The size and luminosity of 
this sample are measured from publicly-available _R-band 
continuum subtracted Ha images using the same methods 
applied to the WiggleZ data as described in Section 3. We 
identified 2091 regions from a total of 7 galaxies. Th e Extra- 
galactic Distance Database (EDD. lTuUv et al.ll2009l ) is used 
for consistent distance indicators for the SINGS galaxies. 

The data for giant HII regions are taken from 
the literature from different parent populations in- 
cludi ng ULIRGS, dwarf irreg ul ars, and isolate d sp 



rals (Gallagher & Hunter! [l983; lArsenauIt fc Rovl 1 19881: 



Rozas et al.ll200d ; iBastian et al1l2006l ; iMonreal-Ibero et all 
20071 '). Studies that included Ha velocity dispersion were 



preferentially chosen from the literature for comparison to 
the high-redshift IFS data. Values are corrected to the dis- 
tances from the EDD and cosmology denned in Section 1. 
Corrections to published data that we introduce to improve 
the comparison of different samples are described in detail 
in the Appendix. 



3.1 Region size 

3.1.1 How to measure size: isophote vs. profile fitting? 

The size of an individual star-forming region is a difficult 
but important physical property to accurately measure. Sub- 
sequent calculations that depend on the size include re- 
gion luminosity, velocity dispersion, relative velocity, and 
metallicity. Two methods are traditio nally employed to mea- 
sure region size: the isophote method (|Sandage fc Tammannl 



Il974 ) and profile fitting, historically called the 'core' method 
( Kennicuttl 19791 ) . In the isophote method a radius is derived 
from the area enclosing a total flux greater than a defined 
fraction of light or surface brightness level. The isophotal re- 
gions selected can be any shape with the radius then char- 
acterised from the total area assuming circular symmetry. 
Depending on the isophote level selected, the regions iden- 
tified can represent either the cores of HII regions or the 
core and surrounding diffuse nebulae. The standard choice 
of isophote varies among studies, with some isophotes se- 
lected to encompass a certain percentage of flux while oth- 
ers are selected visually. Conversely, the core method, or 
profile fitting, measures size by fitting a light profile, most 
commonly a Gaussian, to the surface brightness profile of 
each region. The core radius primarily measures the ionised 
central regions of HII regions, above the surrounding diffuse 
nebulae. 

We measure the isophotal and core radii of HII re- 
gions in seven galaxies from the SINGS survey (NGC0024, 
NGC0628, NGC0925, NGC1566, NGC3938, NGC4254, and 
NGC7552). The isophote method selects regions in the 
SINGS galaxies with surface brightness greater than 1 x 
10~ 16 ergs s" 1 cm -2 arcsec -2 . This isophote was selected 
for the ability to distinguish between close bright HII re- 
gions while still selecting the majority of HII regions in the 
core, arms, and halo of the listed galaxies. Isophotal radii, 
n so , were calculated from isophote areas assuming circular 
symmetry. The core sizes, r coro , were measured by fitting 
a Gaussian to the ID radial surface brightness profile of 
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Figure 1. Comparison between radii and luminosity of HII regions measured from the isophote method and the core method in seven 
galaxies from the SINGS survey (NGC0024, NGC0628, NGC0925, NGC1566, NGC3938, NGC4254, and NGC7552). We find that the 
core method measures systematically smaller HII region sizes (left panel) within a small range of sizes but the luminosity measured in 
both methods (right panel) is consistent. This discrepancy is due to the inclusion of diffuse emission included in some isophotal radii. 
Although diffuse emission can significantly inflate region sizes, region luminosity is dominated by core emission. Dashed black lines show 
one-to-one relations. 



each HII region, which provides a good approximation to 
the data. This profile is chosen because it is the commonly 
used in the literature data, and the signal-to-noise (S/N) of 
the data does not justify a more sophisticated fit. 

Fig- HI shows a comparison between radii measured from 
the isophote method and the core method for the SINGS 
galaxies. Whilst the isophote method measures a larger 
range of sizes (left panel), the luminosity measured in both 
methods (right panel) is approximately consistent. The dis- 
crepancy in measurements arises from how each method ac- 
counts for the Ha background, which can affect both size 
and flux measurements. Small errors, or deviations, in back- 
ground levels produce large errors in isophotal radii, espe- 
cially where the diffuse Ha background is bright. Local dif- 
fuse tails or halos of Ha emission around HII regions can also 
greatly increase the radii causing a significant increase in the 
scatter of measured sizes. Thus, both image quality and H 
II region sub-structures can contribute to the discrepancy 
between sizes from the different measurement methods. De- 
spite the effect of the local background on size, the isophote 
method yields little to no offset to higher luminosities be- 
cause the HII region luminosity is dominated by the core. 
In the core method the total flux of the core is measured 
above the best fit local background level, a free parameter 
in the Gaussian fit, and therefore the core size measurement 
contains minimal bias due to the background emission. 

The issue of background becomes particularly salient 
at high redshift. Diffuse HII halos are either undetected due 
to their low surface brightness or are smeared together in- 
creasing the local background. Therefore, at high redshift 
the local background could become a significant contami- 
nant, increasing luminosity at a given size as a result of 
poorer resolution. This would also be an issue when regions 



are close together, separated by just a few pixels, causing 
clump light to overlap and combine. 

We argue that the core radius is an observationally bet- 
ter determined parameter. The boundary of a core is bet- 
ter defined physically and is less likely to be contaminated 
by the local diffuse nebular gas. No surface brightness level 
needs to be assumed or tuned to select HII regions. This 
method also better serves the comparison of low redshift 
and high redshift data as no flux level or signal-to-noise is 
required for direct comparison. This could be particularly 
important for high-redshift data in which clumps are only 
marginally resolved and often separated by only a few pixels. 
Finally, r COTC is better justified theoretically because it more 
closely approximates the central ionised cluster, which will 
become relevant for future discussions of what drives the re- 
lations between radius and and other observed parameters 
(Section 5). A disadvantage of the core method is that it as- 
sumes the light profile is a 2D Gaussian and will be subject 
to systematic error for other distributions. 



3.1.2 WiggleZ sizes 

Turning now to the WiggleZ sample, we measure the core 
and isophotal radii of 8 clumps from the 3 galaxies which 
show multiple resolved clumps, WK0912.13R, WK1002.61S, 
WK0909_06R. Core radii of the WiggleZ clumps were mea- 
sured by fitting multiple 2D Gaussians simultaneously to the 
unsmoothed Ha images. The measurements given in Table|4] 
represent the 1-sigma Gaussian widths. In cases were clumps 
are elliptical, an average radius was calculated. To measure 
the isophotal radii the isophotes used, 4.5 — 6.5 x 10~ 15 ergs 
s~ cm - arcsec - , were tuned to obtain a reasonable com- 
parison of r COTC and n so . Both methods identified the same 
regions with comparable luminosities. Measurements of r coro 
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Table 4. Properties of Clumps in WigglcZ galaxies from OSIRIS Ha datacubc 



"net 

(km s- 1 ) 



SFR^ 
(M yr" 1 ) 



12 + log{0/R) b 



ID 



u rel 

(km s- 1 ) 



(pc) 



(pc) 



log[M,(M )] 



WK0912.13R 




88.4 ± 4.4 


20.3 ± 6.8 








<8.47 


clump- 1 


55.4 ± 67.4 


140.5 ± 66.5 


2.80 ± 1.87 


1400 


1800 


10.5 


<8.76 


clump- 2 


-39.5 ± 18.1 


52.03 ± 18.1 


1.44 ± 0.64 


1000 


1100 


9.5 


<8.60 


WK1002.61S 




104.3 ± 27.7 


35.6 ± 2.3 








<8.37 


clump- 1 


27.5 ± 6.4 


77.65 ± 6.41 


8.31 ± 0.89 


1600 


1800 


10.0 


<8.42 


clump- 2 


3.2 ± 11.9 


90.74 ± 11.9 


2.54 ± 0.43 


900 


1200 


10.0 


<8.46 


clump-3 


-32.7 ± 22.9 


113.5 ± 22.8 


4.98 ± 1.32 


1700 


700 


10.4 


<8.62 


WK0909.06R 




153.3 ± 8.4 


54.8 ± 3.9 








8.55 ±0.10 


clump- 1 


130.7 ± 7.5 


135.6 ± 7.48 


6.62 ± 0.49 


1200 


2100 


10.4 


8.57 ± 0.10 


clump- 2 


-45.7 ± 6.2 


72.68 ± 6.23 


2.13 ± 0.24 


800 


1500 


0.7 


8.38 ± 0.29 


clump-3 


-88.3 ± 6.4 


80.02 ± 6.43 


11.3 ± 1.2 


3000 


1100 


10.3 


8.44 ± 0.21 



"Velocity relative to the systematic redshift. 
'Calculated within 1 r core . 




Figure 2. Smoothed Ha flux maps of three clumpy galaxies in the WigglcZ kinematic sample observed with the OSIRIS integral field 
spectrograph. Clumps identified in Section 3 are shown by the white ellipses. Pixel size corresponds to 0.05" and the physical orientation 
are denoted by the compasses in the bottom right corner. Images are smoothed with a Gaussian kernel with FWHM = 0.15". 



and ri S o are given in Table [4] and the identified clumps are 
shown in Fig. [2] with the core region sizes overplotted on 
the Ha images. In the remainder of this paper we use the 
core radii of the WiggleZ clumps to represent their sizes, 

^cl -- ^"core ■ 

It is likely that the measured clump radius over- 
estimates the actual size of clump-3 in object WK0909_06R. 
Admittedly this represents a failing of the core method, 
such that low S/N diffuse regions are difficult to model with 
a Gaussian. In object WK0912_13R three clumps can be 
identified visually in the smoothed image (Fig. [2]). However, 
when fitting 2D Gaussians to the raw images the addition 
of a third is not justified by the data as this region is un- 
resolved. This region is also not identified at the selected 
isophote level, which yields similar sizes for clump- 1 and 
clump-2 from the core and isophote methods. Given the un- 
certainty in measuring clump sizes, arising from degeneracies 
in the 2D Gaussian fits, resolution effects, and systematic 
error introduced if the regions are not truly Gaussian, we 
estimate a ~ 30% error on the sizes. 



3.1.3 Resolution effects 

The 8 WiggleZ clumps presented here are resolved as are the 
star-forming regions in the comparison samples. However, 
it is possible that some WiggleZ clumps are aggregates of 
unresolved regions since not all clumps are spherical. The 
resolution of the OSIRIS Hq images studied in the WiggleZ 
sample is ~ 520 — 840 pc in physical coordinates at z = 1.3. 
Because the WiggleZ clumps are close to being point sources, 
the Hq luminosity may be underestimated due to the Strehl 
ratio (~ 30% estimated from the tip-tilt stars) as a result 
of light being spread out by the uncorrected seeing-limited 
halo. 

Despite the regions being resolved, the image resolution 
may still have implications for the scaling relations investi- 
gate in Section 5. Resolution can have the following effects: 
(1) a core and diffuse emission are smeared together, (2) 
multiple core regions are smeared together, (3) diffuse emis- 
sion is smeared to resemble a core, and (4) small regions 
are reduced to noise, with all effects resulting in inflated 
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size and luminosity measurements fo r the r egions observed 
|Pleuss et al J 12000) . ISwinbank et all l^OOO ) investigate the 
effects of redshifting local galaxy images to z > 2 and de- 
grading the resolution on the size and luminosity of regions. 
They find that only the cores of the highest luminosity H 
II regions could be observed at high redshift and that low- 
ering the resolution of images blends multiple regions into 
larger more luminous regions. Both studies find that degrad- 
ing the resolution results in a steepening of the empirically 
measured slope between size and luminosity. We readdress 
these resolution effects in Section 5. 



3.2 Region luminosity 

The luminosity of each WiggieZ clump is measured from 
the integrated spectrum created from the spaxels within the 
1-sigma contour of the best Gaussian 2D fit. A Gaussian 
profile was fitted to the Ha emission in the clump inte- 
grated spectrum with the wavelength, intensity and width 
as free parameters. Errors were calculated by adding in 
quadrature an estimation of the sky noise from the vari- 
ance of the spectrum offset from the emission lines with the 
Poisson error of the photon counting statistics in the emis- 
sion lines. Clump star formation rates (SFR) range from 
1 — 12 Mp yr" 1 derived from the Ha emission using the 



iKennicuttl (J1998T ) conver sion with a Baldry-Glazebr ook ini- 
tial mass function (IMF: lBaldrv fc Glazebrookll200l BG03) 
and are given in Table [4] with star formation rate densities, 
E, = SFR/7T7-;? orc ~ 0.4 - 2 M yr" 1 kpc~ 2 . All the local 
and high-redshift data have been converted to a BG03 IMF 
for comparisons in Section 5. A division by 1.82 is required 
to convert a SFR with a Salpter IMF to a SFR with a BG03 
IMF. 

For WK0912_13R, WK1002_61S, and WK0909.06R the 
clumps make up 21%, 45%, and 37% of the total star for- 
mation respectively with each clump contributing ~ 10%. 
This is in rel ative agreement wi th other IFS studies of 
emission lines (|Genzel et alj|201lf) and broad band studies 
(|Forster Schreiber et al.ll201ll ; lGuo et alJl201lT ). In galaxies 
in the Hubble Deep Field at 1.5 < z < 2.5 up to 50% of 
star forma tion occurs in cl umps with a ~10% contribution 
per clump (|Guo et alj|2011f). In t h e hyd rodyn amic cosmolog- 
ical si mulations of lAgertz et al.l (|2009l ) and ICeverino et al.l 
1)20101 ) the giant clumps account for about half the total star 
formation at a given time and 20% of the disc mass. For com- 
parison, HII regions in H a can contribute 30-50% of the gas 
in lo cal spiral galaxies (jFerguson et al.l 1 1990 ; 



200C ) and up to 75% with beam-smearing 



Pleuss et al 



Pleuss et al 



2000 '). Thus the contribution from clump light measured in 



the WiggieZ sample agrees with high-redshift emission-line 
samples, well-resolved low-redshift samples, and simulated 
samples. 



3.3 Region velocity dispersion 

The velocity dispersion of the gas will be dominated by the 
core region of emission. Although radii measurements are 
critical for size-luminosity comparisons, variations in the 
size of the regions h ave little effect on the integrated ve- 
locity dispersion (e.g. lArsenault fc Rovlll988r ). We test this 
by measuring the velocity dispersion from the integrated 



spectrum in the clumps of the WiggieZ galaxies within 
regions of increasing radii. We find that the velocity dis- 
persion does not change by more than one standard de- 
viation when the radius is increased by a factor of two. 
As such , although some veloc ity dispersions in the local 
sam ple IIArsenault fc Rovll 19881) and the hig h-redshift sam- 
ple (jjones et al.l 12010 ; ISwinbank et al.l 120091 ') are measured 
within isophotal radii, the effect on their velocity dispersion 
measurements should be minimal. 

Velocity dispersions, <T nct , for the WiggieZ clumps are 
derived from the integrated clump spectrum and corrected 
for instrumental broadening by subtracting the OSIRIS res- 
olution in quadrature and are given in Table [4] We use this 
measure of velocity dispersio n rather than the f lux-weighted 
veloc ity dispersion, cr mC an (ILaw et al.l 120071; iGreen et al.l 
120101 ) or other measurements ( Davies et al.ll2011f ). for bet- 
ter comparison to the values taken from the literature as it 
more closely resembles a long-slit observation. Rotation sig- 
natures across the clumps are measured to be ~ 20 — 80 km 
s _1 , with average Av/Ar c i ~ 15 km s" 1 kpc -1 , where At> 
is the change in velocity across the clump. This is a small 
effect in comparison with measured velocity dispersions, as 
such no correction has been applied. Velocity dispersions 
of the clumps, uncorrected for beam-smearing, are either 
consistent with or slightly lower than the total disc disper- 
sions. Total disc dispersions, measured in Paper I, are given 
in Table [3] Recent simulations show that even with mod- 
erate beam smearing — 0.1 arcsec — the actual observed 
clump dispersion is highly contaminated by the disc disper- 
sion making it difficult to d istinguish between the two signa- 
tures (|Ceverino et alJl201lT ). However, the same simulations 
predict that clump and disc dispersions should be physically 
comparable in the absence of beam- smearing, minimising 
the ramifications for our measurements. 

Clumps that have higher mass surface densities than the 
disc, E c i >> Ed, could lead to higher dispersions within the 
clumps due to intense star formation or stellar winds. Large- 
scale winds can be launched from regions forming stars at 
E* > 0.1 M y^ 1 kpc" 2 (JHeckman et all 120001 ) and in- 
deed the WiggieZ clumps measured here have star formation 
surface densities > 10 x this threshold. If large-scale stel- 
lar outflows are driven by the star formation in clumps, in- 
flated velocity dispersion s are expected due to an underlying 
broad wind component ([Heckman et al.lll990 ; iGenzel et al.l 
120111 ). In Fig. [3] we stack the integrated spectrum from all 
8 WiggieZ clu mps at rest waveleng th, removing any veloc- 
ity broadening (jShapiro et al.ll2009l ). The stacked spectrum 
is fit by two Gaussian components, a narrow peak (FWHM 
~ 170 km s" 1 ), and underlying broad peak (FWHM ~ 490 
km s" 1 ) which has a consistent systemic velocity with rela- 
tive offset —24 ± 30 km s -1 . We find that a 2-Gaussian fit 
provides a significant improvement in the value of the chi- 
squarcd statistic compared to a 1-Gaussian fit, producing a 
decrease Ax 2 ~ 60 for the addition of 3 extra parameters 
and a reduced \ 2 °f 1-07. The broad Gaussian contributes 
40% to the total flux. We rule out AGN as the driver of 
the broad Ha component because the clumps are off-centre 
from the photometric and kinematic centre of the galaxies 
(jAlexander et alj|2010n and no evidence is seen on galactic 
scale for AGN contamination, as discussed in Paper I. 

The broad velocity component observed in the stacked 
WiggieZ clumps has a smaller velocity width and smaller 
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ing [PettinrTPagel |2004ft ; 
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Figure 3. Stack of the integrated Ha spectra from 8 WigglcZ 
clumps corrected for velocity broadening. The red line is the best- 
fit model to the Ho and [Nil] emission lines composed of two 
Ha Gaussian components. The blue dotted line is the narrow 
component with FWHM = 170 km s — 1 and the blue dashed line 
is the broad component with FWHM = 490 km s _1 which has 
a consistent systemic velocity with relative offset — 24 ± 30 km 
s _1 . The broad component contains 40% of the total flux. The 
orange line is the residual of the best-fit model subtracted from 
the data. The narrow black dotted line at zero velocity shows the 
instrumental resolution from a stack of sky lines in the observed 
spectra. Vertical dotted lines show the expected positions of [N 
IIJA6548, HaA6563, and [NII]A6584. 



blueshift compared to z ~ 2 starbursts galaxies, which may 
have broad components wit h up to FWHM ~ 1500 km s -1 
blues hifted by 50 km s^ 1 (JGenzel et al.ll201ll : Friran et all 
I2011T ). However, we find that the broad velocity compo- 
nent of the WiggleZ clumps contributes a comparable or 
larger fraction to the overall flux, equal to 40%. Although a 
broad component is not detected in the spectra of individual 
clumps, the presence of a significant broad component in the 
stack of clump spectra indicates that individual clump ve- 
locity dispersion measurements may be over-estimated due 
to contamination from large-scale winds. By comparing to 
a single Gaussian fit to the Ha emission of the clump stack, 
this is estimated to be approximately a 15% effect, or on 
average ~ 14 km s _1 , comparable to the errors given in Ta- 
bled The possible effect of winds on both clump luminosity 
and velocity dispersion is reflected in the error bars in the 
subsequent figures. 



3.4 Clump metallicity and ages 

It has been shown that the scatter in HII region mea- 
sure ments can be partially explained by region m etallic- 
ity (JTerlevich fc Melnickl flfjJH ; iMelnick et ail 1 1987ft . Vari- 
ations in metallicity may also reveal details about the age 
of clumps, important to their evolution. We calculate clump 
metallicities for each individual WiggleZ clump and for the 



(1) 



12 + log f g J = 8.90 + 0.57 x log f H3 j . 

These values are given in Table 2] For individual WiggleZ 
clumps a single Gaussian is fit to the Hq and [Nil] emission 
lines. For the clump stack, the Ha flux used in equation ([TJ is 
from the narrow velocity component (shown in Fig. [3] as the 
dotted blue line). The metallicity of all but the WK0909.06R 
clumps are upper limits as [Nil] in this sample is typically 
undetected. The stack of eight clumps has a metallicity of 
12 + log(0/H) = 8.35 ± 0.31. All metallicities in individual 
galaxies are consi stent within the erro rs, w hich agrees with 
the sim ulations of llmmeli et al.l (|2004aft and lBournaud et al.l 
(J20081 ) and implies a common clump evolution history. 



4 DISC INSTABILITY 

As outlined in Section 1, models predict that clumps form 
out of gravitational instabilities in gas-rich turbulent discs. 
The high levels of turbulence are predicted to drive up the 
expected instability scale in disc galaxies resulting in mas- 
sive star-forming regions. In this section we characterise the 
size and mass instability scale of the WiggleZ discs given 
the high velocity dispersions measured in Paper I and the 
previous section. We compare the WiggleZ clumps to the 
high-redshift and low-redshift samples, which have velocity 
dispersions ranging from 10— 100 km s . We calculate rela- 
tive Toomre parameters, shown in Fig. [4j to identify regions 
of instability where clumps could form. 

For star-forming regions to form in gaseous discs, the 
galaxies must become gravitationally unst able such that 
Q < 1 where Q is the Toomre parameter IjToomrel ll964T ) 
defined for a pure gas disc as 



Qe 



7rGE e 



(2) 



To test these models of instability we calculate Q gas for the 
three galaxies hosting clumps in the WiggleZ sample using 



Qi 



Vc /g 



(3) 



as derived from lGenzel et al.l (|201lh . where k is the epicyclic 
frequency, a is the local disc gas velocity dispersion, v c is the 
circular velocity, and E g is the local gas surface density in 
the disc. The constant a relates the epicyclic and circular 
frequency in different poten tials such that a = x/2 for a disc 
with a flat rotation curve (JBinnev fc Tremainj 120081 ). We 
estimate the gas fraction, f g , from the inverse Kennic utt- 
Schmidt law fKS: ISchmidlll95filKennicutt et al.ll2007r ). 

In Paper I flat disc models were fit to the Ha velocity 
maps of the WiggleZ galaxies. We obtained estimates for 
v c sin(i) rather than v c because the inclination could not be 
constrained for lack of deep broadband data and the full ex- 
tent of the kinematic gas disc is not observed. As a result 
the actual value of the Toomre parameter across the disc 
is highly uncertain. A relative Toomre parameter however 
is sufficient to identify the relative stability of the galaxy 
discs. We note that within the errors of v c sin(i) for the 
disc models, Q sas for each galaxy has regions of instabil- 
ity (Qgas < 1). In Fig. [4] we show relative Toomre maps 
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Figure 4. Toomrc parameter, Qgas (equation [3} , mapped over three clumpy galaxies in the WigglcZ kinematic sample. The calculation 
of the Toomre parameter assumes the Kennicutt-Schmidt law and is uncertain due to the degeneracies of the variables in equation J3j). 
As a result a relative linear scale for Q gas is shown. All three WigglcZ galaxies have areas that are consistent with being in an unstable 
state (Q gas < 1) within the errors. Clumps identified in Section 3 are shown by the white ellipses. Likely unstable regions are shown in 
purple and dark blue. The clumps correlate approximately with Toomre unstable regions. 



for WK0912.13R, WK1001.61S, and WK0909.06R. Over- 
laid on the maps are the clumps identified in Section 3 from 
the core method. The clumps identified from the Ha maps 
match to the less stable regions of the discs as expected 
from the models. These results are in agreement with IFS 
data at z ~ 2 i n which Toomre in stability correlates with 
clump location (JGenzel et al.ll201lr ). It is expected that low 
values of Q gas will correlate with regions with high SFRs 
as the Toomre parameter is proportional to the inverse of 
star formation as a result of how the gas fraction is calcu- 
lated. However, significant peaks in local velocity dispersion 
coincident with the clumps (for example if the clumps are 
virialised and rotationally supported) could counterbalance 
the local reduction of Q gas by the star formation. However, 
this effect is not observed since significantly higher velocity 
dispersions are not measured in the clumps versus the disc. 
When discs are gravitationally unstable, clumps are ex- 
pected to form under Jeans collapse on a characteristic scale. 
For the WiggleZ galaxies the expected masses and sizes of 
the clumps would be limited by the Jeans mass and length, 
given by 



and 

A/j = 



m 



4tt /Xj 

t It 



1/2 



(4) 



(5) 



where <Td is the average velocity dispersion in the discs es- 
timated by the velocity dispersions of the clumps and H 
II regions, and p is the density of the disc. By combin- 
ing equations ((4]) and (0) and assuming that the mea- 
sured clump diameter is equal to a Jeans radius, we cal- 
culate Jeans masses for the WiggleZ clumps of Mj = 
1.1 x 10 9 — 1.1 x 10 10 M fl , which are in agreement with the 
models (e.g. lKrumholz fc Dekel2010r ). In Fig.[5]we compare 
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Figure 5. Ha Luminosity vs. Jeans mass for the WiggleZ clumps 
(red diamonds), high-redshift comparison sample (blue filled sym- 
bols) and low-redshift giant HII region sample (gray points). 
Ha luminosity and Jeans mass (Mj = 7r 2 ir 2 r c i/6G) show a 
tight correlation for star-forming regions in local star-forming 
galaxies (HII regions) to turbulent high-redshift unstable discs 
(clumps) with sizes of ~ 50 — 5000 pc. The blue dashed line is 
the least squares best fit to the data, given by Z/jj a [ergs s_1 ] = 
8 x 10 29 Mj[M ] 124±0 - 05 . The black dot-dash line represents 
Ln a oc Mj arbitrarily normalised for comparison. 

the Jeans mass to Ha luminosity for the WiggleZ clumps 
(red diamonds), the higher-redshift clumps (blue squares 
and triangles), and the low-redshift giant HII regions (gray 
points). A tight continuous correlation is seen for all star- 
forming regions from Mj = 10 6 — 10 10 Mq. The correlation 
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Figure 6. Relations between the Ha luminosity, size, and velocity dispersion of HII regions in local star-forming galaxies and clumps 
in high-redshift galaxies. Clumps fro m WiggleZ galax i es are represented by the b right red diamonds. High-redshift clumps are shown as 
filled symbols in blue and taken from lSwinbank et al.l J2009h ; [jones et all }201ol) ; lGenzel et all l Effllil) ; iForster Schreiber et alj (l201lh . H 
II regions from loc al SINGS galaxies ar e displ a yed as black dots, giant HII regions from loca l galaxies are shown as open black points 
and a re taken from lGallaeher fc Hunterl jl983ft ; lArsenault fc Roy| jl988h ; lBastian et all l l2006l) ; lRo"zas et all l l2006l) ; lMonreal-Ibero et all 
(2007). The blue dashed lines show the least squares best fits, given in the top left corner of each panel with correlation coefficients, R , 
given directly below. Simulated clumps at 1.9 < z < 3.0 from lCeverino et all KOllh are plotted in orange but have not been included in 
the fit. 



in Fig. [5] combines the three main observables discussed in 
this paper with scatter characterised by a correlation coeffi- 
cient of 0.90 and least squares best fit of Lu a oc A^fj ' 24 , 
indicating that the mass closely traces the ionised gas of 
star-forming regions. 

5 SCALING RELATIONS OF STAR-FORMING 
REGIONS 

Are the kpc-sized clumps found at high redshift formed out 
of disc instabilities in a different manner to local HII re- 
gions, or are they more massive analogues as suggested by 
Fig. [Sf We show the correlations for measured Ha luminos- 
ity, velocity dispersion, and size for the WiggleZ clumps, the 
high-redshift comparison sample, and the low-redshift com- 



parison sample in Fig. [6] We find tight correlations between 
these properties over 50-5000 pc for HII regions and clumps 
taken from eleven different studies. Dispersions range be- 
tween ~ 10 — 100 km s" 1 and Ha luminosities span five or- 
ders of magnitude. Assuming equal weighting for all points 
we use linear least squares fitting to find the best fit to the 
relations to be 



log(o-) = (0.42 ± 0.03) x log(d) + (0.33 ± 0.09), 
log(L HQ ) = (2.72 ± 0.04) x log(d) + (31.99 ± 0.08), 
log(L HQ ) = (4.18 ± 0.21) x log(o-) + (33.61 ± 0.31), 



(6) 



where Lua is measured in ergs s~ , d in pc, and the a in km 
s _1 , with correlation coefficients, R 2 , of 0.72, 0.82, and 0.87 
respectively. These fits are overplotted as blue dashed lines 
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in Fig. [6] Simulated clumps are included in Fig. [6] (orange 
points) for comparison to the data, but are not included in 
the above fits. Below we discuss the scaling relations and 
their possible drivers. 



5.1 Velocity Dispersion vs. Size 

We derive the relation between size and velocity dispersion 
for regions that form out of Jeans collapse in an isothermal 
disc. The disc scale height, H, for an isothermal disc is given 
by 



H=- p(z) dz, (7) 

Po Jo 

(|Spitzerlll942l ) where po is the midplane density value 
and z is the height above the plane. The scale length has 
the solution for a purely self-gravitating unmagnetised gas 
dominated disc 



H 



(2tyG Po o)^ 2 7rGE g ' V ; 

where H^ and p^ are the disc scale height and den- 
sity in this limit with a defined as the turbulent veloc- 
ity dispersion. The surface density of the disc is given by 
Eg = 2H 0O p O2 . Therefore we obtain an equation for the den- 
sity such that 



Poo 



ttGS 2 
~2J r 



(9) 



Combining this with the definition of the Jeans length 
in equation (|4]), assuming p = poo, we obtain a threshold 
length for the regions forming under Jeans collapse in an 
isothermal disc such that 
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Figure 7. Star formation surface density (£*) vs. gas mass sur- 
face density (S g ; derived from equation HOP for the same data 
described in Fig. [6] with all giant HII regions now shown by gray 
points. The bla ck dot-dash line is the Kennicutt-Schmidt law, 
£* oc Si 4 from iKennicutt et al.l IJ2007I ). The average gas surface 
density for the WiggleZ clumps, assuming the KS law, derived 
from the measured S* is ~ 500 Mq pc -2 in comparison to what is 
derived from the Jeans equations for a pure gas disc of E g ~ 1000 
Mq pc~ 2 . The blue dashed line is the least squares fit to the data 
given by £»[M yr" 1 kpc~ 2 ] = 2 X 10~ 5 S g [M pc" 2 ] 1 ' 37 * - 14 . 
The horizontal line represents the threshold of star formation 
required to host stellar winds IIHeckman et al.l 120001 ), The star 
formation surface density and KS law are plotted for a Baldry- 
Glazcbrook IMF. 



Aj< 



GE E ' 



(10) 



The corresponding mass limit from equation (JSJ) be- 
comes M] < v / 2i" 2 o" 4 /6G 2 E g . To associa te this scaling with 
measured size we adopt the convention of lDekel et al.l (|2009l ) 
which relates the typical clump radius to the Jeans scale as 
r c i ~ Aj/2. Thus, we define region diameters to be ~ Aj as 
plotted in Fig. [6] We fit equation (|10|l to the first panel in 
Fig. [6] and obtain a normalisation for the gas surface density 
of E g =6OOM pc" 2 . 

We use the correlation between size and dispersion de- 
fined in equation (|10[1 to estimate the gas surface density and 
test the Kennicutt-Schmidt law for individual star-forming 
regions. We assume that at the time of collapse, star- forming 
regions should have comparable surface densities and veloc- 
ity dispersions to the disc in which they formed, and that 
the mass of the star-forming regions is dominated by gas 
(Eg = v2cTci/Gr c i)- This is only an approximation as the re- 
gions are not purely molecular. Therefore, it is expected that 
the gas surface density determined by equation (jlOp should 
always be higher than the molecular gas surface density as 
estimated from the KS law. In Fig. [7] the star formation 
surface density is plotted against the derived gas surface 
density for the low an d high-redshift regions. The KS law 
(jKennicutt et al.ll2007l ) is overplotted by the black dot-dash 



line corrected to a BG03 IMF. The derived gas surface den- 
sity for all regions is indeed greater than expected from the 
KS law, but is otherwise a reasonable approximation to the 
true gas surface density. The linear least squares fit to the 
data gives a slope of 1.38±0.14 as shown by the blue dashed 
line and correlation coefficient of 0.66, consistent with the 
slope of 1.4 predicted by the KS law. We note that the fit 
degrades at high redshift. This could be to a selection effect 
given that we can only observe regions with the highest star 
formation surface densities. 

Using this method we obtain upper limits on the gas 
surface densities of order ~ 500 — 3000 Mq pc~ 2 for Wig- 
gleZ clumps, which overlap with the star formation surface 
densities for giant HII regions (~ 100 — 3000 M pc -2 ), and 
are comp arable to the mean gas surface density fo r clumps 
at z ~ 2 (jCeverino et al.ll2010l : iGenzel et al.ll201ll ). To em- 
phasise the difference between the local and high-redshift en- 
vironments in these instances, clumps have been compared 
to local giant molecular clouds (GMCs) w hich have typica l 
surface densities of ~ 60 - 170 M pc ~ 2 (JBlitz et al.ll2007i ; 
iBolatto et al.l 120081 : Irlever et a.1.1 l2009h . However, we argue 
that clumps are more comparable observationally and phys- 
ically to HII regions than GMCs as they are observed from 
their OB stars, which have already formed out of molecular 
gas. 

The size-velocity dispersion relation shows the most 
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scatter in Fig. [6] with correlation coefficient of 0.72. Is this 
scatter due to measurement errors or an expected spread 
in region surface density in accordance with the KS law? 
The scatter in velocity dispersion ranges from ~ 20 — 90 km 
s _1 while the average error on the measured velocity dis- 
persion of the WiggleZ clumps is ~ 20 km s _1 . This scaling 
is expected to show scatter due to measurement errors dis- 
cussed in Section 3: size is the most uncertain parameter 
to measure, and the contribution of galactic winds in and 
around star-forming regions may result in added scatter in 
the velocity dispersion measurement. However, if regions are 
expected to follow the KS law (as Fig.[7]tentatively suggests) 
then from the three order-of magnitude spread measured in 
E* a range of ~ 4000 Mope -2 in E g would be expected. 
This range translates into a spread in velocity dispersion of 
~ 60 km s at a given radius, approximately equivalent to 
the average scatter in Fig. [6] and greater than the typical 
errors on velocity dispersion of the WiggleZ sample. Thus, 
if the KS law holds, then the expected variations in the gas 
surface density can explain the scattered relation between 
size and dispersion. 

The same scaling between size and velocity dispersion, 
a oc r 1 , can be obtained for virialised objects from the 
equation of dynamical mass (Md yn = 5<r 2 r c i/G) and assum- 
ing a constant surface density such that a 2 = 7rGEr c i/5. 
When this relation is fit to the data in Fig. [6] the best-fit 
average mass density for all regions, E = 600 M s pc~ 2 , 
is equivalent to the surface density derived using the Jeans 
scaling. From the scatter in Fig. [6] the mass densities for all 
regions ranges from 100 — 5000 Mq pc -2 similar to what 
is seen in Fig. and discussed above. If regions have viri- 
alised then their mass will be better represented by the dy- 
namical mass, which follows the same scaling as the Jeans 
mass (M oc a 2 r/G) determined in Section 4. However, the 
virial scaling yields masses a factor of 3x greater than their 
Jeans mass. When the dynamical clump masses are added 
within their respective discs they approximate ~ 50 — 90% 
of the total stellar mass, compared to the Jeans masses 
which amount to ~ 20 — 30% of the total stellar masses 
of the discs. The Jeans masses are more comparable to the 
expected masses from individual clumps, and to the total 
mass contributi on of clumps to the whole system, predicted 
by simulations JDekel et al.ll2009l ; IKrumholz fc DekelllJOlol ; 
ICeverino et al.l 120111 ). Thus, it seems unlikely that all the 
observed clumps have virialised and constitute such a large 
mass fraction of their discs. 

We note that the high-redshift clumps, from the Wig- 
gleZ and SINS datasets, on average have higher velocity dis- 
persions than predicted by the best-fit relation described by 
equation (6). If we assume that the clumps from both sam- 
ples are con taminated by galact ic scale winds, as suggested 
here and in lGenzel et al.l (|2011f ). with an average contribu- 
tion comparable to the result of our stacking analysis in 
Section 13.31 then the high-redshift dispersions are reduced 
by a factor of ~ 0.85 making them more consistent with the 
scaling relation of equation (6). 

5.2 Luminosity vs. Size 

The scaling of Ha luminosity with size of HII regions is 
well understood if size is represented by the radius of a 
Stromgren sphere, r s , the boundary between ionised and 



neutral hydrogen. The luminosity of an ionisation-bound re- 
gion scales as Lhq oc r 3 , or more specifically 



4t N H 2 a B /ic 



A, 



(11) 



where Nh is the ionised hydrogen number density and qb 
is the Ca se-B recombinatio n coefficient assuming Tun = 
10 000 K l|Osterbrockill989l ). We obtain a density of N H ~ 
500 cm -3 by normalising the scaling relation to the data and 
assuming r B = r c \. For comparison, Nh ~ 100 cm -3 in lo- 
cal ULIRGS (|Monreal-Ibero et alj|20u7h . The residual from 
subtracting a least squares fit, fixing L oc r 3 , from the data is 
shown in Fig. [5J However, this is an idealised relation which 
assumes direct observations of individual Stromgren spheres 
at a constant number density and temperature. The empir- 
ical relation observed in our data yields a slightly shallower 
slope of 2.72 ± 0.04, which could result from a variety of 
factors that act alone or in aggregate. Below we investigate 
the physical and observational processes that could affect 
the slope of the relation as well as consider the scalings of 
GMCs. 

If the hydrogen atoms within the cloud can recombine 
faster than they can be ionised then not all the ionised pho- 
tons will be absorbed. In this case, the ionisation region is 
said to be density bound, as the observed luminosity is lim- 
ited by the density of the cloud, and could explain deviations 
from the scaling relation such that regions are less lumi- 
nous at a given r adius. A detailed discussion of this scenario 
can be found in iBeckman et al.l (|2000f ) and the references 
therein. In this study all regions with Lu a > 10 38 ' 6 ergs s _1 
are predicted to be density bound due to a constant surface 
density for GMCs from which they form. We note that all of 
the WiggleZ clumps and the majority of giant HII regions 
studied here are more luminous than this threshold, which 
suggests that a shallower slope could be observed as a result. 
An important consequence of regions being density bound is 
that escaping photons could be responsible for ionising the 
diffuse interstellar medium in the disc. 

As discussed in Section 3.1.3, resolution can affect the 
size-luminosity scaling relation. Beam-smearing results in 
a steeper slope i n size-luminosity space than expected fo r 
resolved regions (|Pleuss et alJl200Ct ISwinbank et aHl2009h . 
Therefore, if the Stromgren sphere scaling is the correct 
theoretical prediction for the empirical relation of size- 
luminosity, then resolution is not the dominant effect on the 
slope and may counterbalance other processes that result 
in a shallower slope (e.g. when regions are density bound). 
Dust, metallicity, and magnetic fields are also likely to re- 
duce the observed ionised radiation at a given size, however 
a detailed discussion of their effects is beyond the scope of 
this paper. 

If the more appropriate analogues to Ha clumps are 
GMC s rather than HII regions as suggested in th e litera- 
ture (|Murravet alj|20ld : IKrumholz fc ' Dekdl20ld ) then a 
different scaling is expected between luminosity and size. 
For this case, a molecular tracer of radiation is necessary. 
The scaling between CO luminosity and size for molecu- 
lar clouds follows a shallower relation, than expected from 
Stromgren spher es, of order L oc r 5 ' 2 (|Solomon et al.lll987l ; 
lElmegreenl ll989T ). However, although this scaling provides 
a reasonable fit to the HII region data shown in Fig [HI it 
assumes a CO luminosity and cannot be tested by our Ha 
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Figure 8. Residual luminosity vs. velocity dispersion and diameter from linear least squares fits assuming the scaling relations between 
luminosity and velocity dispersion and luminosity and size denoted in the upper left corner of each panel. Symbols and colours are 
the same as in Fig. [6] The relations arc derived for (top) ionisation bound Stromgren spheres that formed from Jeans collapse in an 
isothermal disc, (middle) GMCs, and (bottom) virialised regions. The L — cr relation is best approximated by L oc <r 4 and the L — r 
relation is best approximated by Locr 3 . These scaling relations are discussed in detail in Sections 5.2 and 5.3. The most deviant points 
in the right panels, offset by two orders of magnitude in luminosity, are from locally LIRG merger of NGC4038 and NGC4039 (open 
inverted triangles), lensed sub-L* galaxies at z ~ 2 (blue closed circles), and an L* galaxy at z ~ 5 (blue closed squares). The offset 
could be due to a difference in how regions are measured, as both high-redshift studies used the isophote method rather than profile 
fitting and the determination of sizes in the low-redshift sample is unclear. 



data. We note that a similar relationship, of order L oc r' ! , 
has been observed between lensed clumps at 2 = 2.3, GMC 
core s, and young HII regio ns for L26O , a tracer of molecular 
gas (jSwinbank et al.ll2010f ). These results suggest that the 
underlying physics of the star-forming processes is similar 
from z — 2.3 to the present-day Universe. 

Some observations (jSwinbank et al.l 12003 ; Ijones et al.l 

[201fj) of high-redshift clumps show deviations from the em- 
pirical luminosity-size relation. In these results the clumps in 
both a lensed L, galaxy at z ~ 5 and lensed sub-L* galaxies 
at z ~ 2 are offset by two orders of magnitude in luminos- 
ity from local HII regions (blue filled circles and squares in 
Fig. [6]and Fig. [SJ. An interpretation of the observed offset is 
that these regions, as well as giant HII regions, form due to a 
starburst mode or due to lower metallicities of high-redshift 
galaxies. In comparison, there is a half order of magnitude 
offset in luminosity of the WiggleZ and SINS clumps from 
the best-fit relation to all the regions. A fraction of the offset 
in the WiggleZ and SINs clumps is likely due to an inflation 
of the clump luminosities from large-scale winds. The anal- 
ysis in Section [3~3l suggests that winds can contribute up to 
40% of the luminosity of a region. If the luminosity of the 
WiggleZ and SINS clumps are corrected for this effect the 



remaining offset is less than half an order of magnitude from 
the best-fit relation. While a slight offset in luminosity exists 
for the WiggleZ clumps, it is small in comparison to the two 
order of magnitude offset found for the lensed clumps. The 
remaining offset could be a result of the factors discussed 
in this section. We find no significant luminosity offset be- 
tween HII regions and giant HII regions at a given size. 
Furthermore, we include data from other low-metallicity re- 
gions with comparable values and find that they fall on the 
main relation. Thus, another basis is required to explain the 
deviation of the lensed regions. One possible explanation for 
the discrepancy is that the sizes of the regions from the two 
lensed studies are measured by the isophote method which 
could overestimate the luminosity for a given region by in- 
cluding surrounding flux of the disc. 



5.3 Luminosity vs. Velocity Dispersion 

Assuming from our discussion above that HII regions and 
clumps are idealised Stromgren spheres (L oc r 3 ), which 
form out of Jeans collapse in an isothermal disc (r oc a 2 ), 
then it is expected for luminosity to scale with dispersion 
as L oc a 6 . This scaling can also be derived from the gas 
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accretion rate as a function of viri al velocity and the energy 
deposition rate due to accretion (|Dekel fc Birnboiml 120081 ; 
iLe Tiran et al.ll201ll ). However, as outlined in the previous 
section and shown in Fig. [8] the regions observed are not 
idealised Stromgren spheres and therefore a shallower slope 
is expected for the luminosity-velocity dispersion relation. 
An un- weighted fit to the data reveals log(L) = (4.18 ± 
0.21) x log(cr) + (33.61 ±0.31). 

We re-address the argument that the regions are viri- 
alised, given that the non-weighted fit most closely approxi- 
mates L oc a 4 as the best description of the data. This scal- 
ing is expected from the virial theorem assuming a constant 
mass to light ratio (M IL) and constant surface brightness 
(jFaber fc Jackson! Il976l ). However, as M/L is not constant 
(as shown in Fig [5]) deviations from the predicted scaling 
relations are expected and L oc a 4, is no longer an appro- 
priate test of virialisation. Furthermore, the same scaling, 
assuming M/L is constant, can be derived from the Jeans 
equations such that Mj oc o- 4 /G 2 E g . 

Finally, a relationship between luminosity and velocity 
dispersion of galaxies and th e interstellar medium has re 
cently been observed local ly (IDib et al.l 20061; iGreen et al 



201C ) and at high-redshif t (jLehnert et al 



2009: 



Green et al 



20101 ; iGenzel et al.l 120111 ) . These observations have led the 



authors to conclud e that star formation drives the large ob- 
served line- widths. iLehnert et al.l (|2009P l suggest that the ve- 
locity dispersion is due to mechanical energy released by the 
star formation in the form a — (eE*) ' 5 where e is the effi- 
ciency of star formation and both a and E* can be global 
or local properties. The relationship between star formation 
surface density and velocity dispersion for the star-forming 
regions studied here is shown in Fig. [5] with the best linear 
least squares fit of a - (100±70 x E„)°' 6±0 - 1 represented by 
th e blue dashed l ine. T he fit is consistent with the results 
of lLehnert et al.l (|2009h . However, the fit is also consistent 
with the expected star formation density in a marginally 
stable galactic disc of constant circular velocity assuming 
the KS law on gala c tic an d sub-galactic scales, as noted by 
iKrumholz fc Dekell (|2010r i. In this scenario, placing the KS 
law in equation ((2]), assuming Q — 1, yields a oc 7tGE2' 7 /k. 
This interpretation is more consistent with our results in 
Section 5.1 and assumes the KS law holds for star-forming 
regions (Fig. [7]). However, we note that the a — E* relation 
shows a much greater scatter (_R 2 =0.57) than the relation 
between luminosity and velocity dispersion (_R 2 =0.87). 

In conclusion, we find continuous scaling relations at 
high and low redshift for star-forming regions between Ha 
size, luminosity, and velocity dispersion. When the three ob- 
servables are combined by comparing Ha luminosity with 
the Jeans mass (M oc cr 2 r c i/G) the scatter is significantly 



B? = 0.87, L 



R 2 = 0.82, and 



reduced (L — a 
L — Mj : R 2 = 0.9). Some of the scatter observed in the size- 
velocity dispersion relation may be a natural consequence 
of the Kennicutt- Schmidt law such that a 2 oc E g . The data 
supports the theoretical arguments that star-forming regions 
form out of Jeans collapse in an isothermal disc with their lu- 
minosities approximated by Stromgren spheres, with clumps 
accounting for 20-30% of the stellar mass of the discs. These 
results indicate that star-forming regions at high-redshift 
form through similar processes as star-forming regions lo- 
cally. 
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Figure 9. Star formation surface density vs. velocity dispersion 
for the same data described in Fig. \E\ with all giant HII re- 
gions shown by gray points. The blue dashed line is the least 
squares best fit to the data, given by <r[km s _1 ] = (100 ± 70 X 
S,[Mg yr — 4 kpc - 2 ])0.6±o.i Although this relationship is more 
scattered than the Ljj a — cr relationship shown in Fig.[6]a correla- 
tion is seen, which could be a result of mechanical energy r eleased 
by the star formation (cr = feSt) 0,5 ; ILehnert et al.ll2009T ) or the 
expected star formation density in a marginally stable galactic 
disc of constant circul ar velocity assuming the K ennicutt-Schmidt 
law (<r oc 7tGS9- 7 /k; IKrumholz fc Dekell I20T0T) . The star forma- 
tion surface density is plotted for a Baldry-Glazcbrook IMF. 



6 CLUMP STABILITY AND SURVIVAL 

It may be predicted as a consequence of the changes in 
global environment that star-forming regions at high redshift 
should not resemble star-forming regions locally. The high 
accretion rates predicted to drive up turbulence in the early 
universe are expected to dissipate on cosmological timescales 
for the most mass ive systems and contribu te minimally at 
the present epoch (|Dekel fc Birnboimll2006l h Yet despite the 
increased sizes and dispersions, the regions presented here 
(z ~ 1.3) have overlapping star formation surface densi- 
ties and mass densities with giant HII regions in the local 
universe and form continuous scaling relations with local 
star-formation regions. This indicates that while the turbu- 
lence of discs may have important implications for the size 
and luminosity of regions which form within them, the same 
processes govern their formation from high redshift to the 
current epoch. 

Given the supersonic line-widths observed in the Wig- 
gleZ clumps and local HII regions, star formation is ex- 
pected to occur rapidly, depleting the available gas for star 
formation within a dynamical time. For the gas to remain 
turbulent as is observed, a m e chanism is required to stir 
up the ISM (|Stone et al.lll998l ; [Mac Low et al.lll998l 'l. Yet, 
there is no consensus in the field as to what generates gas 
velocity dispersions across these different epochs. The pos- 
sible drivers of the turbulence at low-reds hift are gravita- 
tional support (|Terlevich fc Melnickl Il98l[ ) , hydrodynami- 
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cal fl o ws in the form of stella r winds (| Gallagher fc Hunterl 



1 19831 ; lArsenault fc Rovl Il988l; iTenorio-Tagle et all Il993h . 
champagn e flows (IGallagher fc Hunterl 1 19831'). supernov a 
feedback JMac Low fc Klessenl |2004 Jjoung et alj l2009h . 
and bow shocks 1 Tenorio-Tagle et al.l 1 19931 ). These gen- 
eral scenarios have also been put forward to account 
for what maintains and drives up the turbulence at 



orted locally by a study of lumino us galaxies at z ~ 0.1 
Green et al. | 120101 ) and in models (jElmegreen fc Burkertl 



high r edshift as well (llmmeli et al.l 2004bl; Bournaud et al.l 



20071; iForster Schreiber et al.1 120061; iGenzel et al.1 1200*1 



Dekel et alj 120091; iKrumholz fc Dekell l20ld; lAumer et ail 



20ld ; lElmegreenfc Burkertl l20ld; iGenel et all l20ld) "in 



addition to gas accretion (|Forster Schreiber et al. 20061: 



Genzel et ai1l200d ; lElmegreen fc Burkertll2010l ; IGenel et al 
201C) and star formation ( Lehnert et al.1 120091 ; iGreen et al 
2010J). In this section, the available data and models are 



compared where possible to explore the implications of the 
previous sections in the context of clump formation and life- 
times. 

Models predict that the dominant feedback mechanism 
at high redshift is likely radiation press ure which can drive 
up the turbulent motion s in the ISM (jKrumholz fc Dekell 
120101 ; iMurrav et al.ll20ld ). Radiation feedback can provide 
pressure support to the clumps as well as drive large stel- 
lar winds. However, if strong enough, radiation feedback 
can act to disrupt the clumps before virial ising and mi- 
grating to the centre to form a stellar bulge II Murray et al.l 
I2OIO1 ; IGenel et al.ll20ld ). Injection of energy into the inter- 
stellar medium by radiation pressure could be confirmed 
by detecting massive radia tion-driven outflows from clumps 
IjKrumholz fc Dekell |2010|). In Section 13.31 an underlying 
broad component is detected, which contributes ~ 40% to 
the total emission line flux from the clumps. The broad 
Hq component represents feedback due to thermal pres- 
sure from ionised gas and shocked stellar winds. These re- 
sults, in combination with similar findings from 2 ~ 2 
clumps (jGenzel et al.1 120111 ) , provide evidence that winds 
from clumps could inject energy into the system, which 
would in turn drive up the turbulent motions in the ISM 
in accordance with the theoretical models. 

If the clumps do survive feedback and clump-clump 
interactions then they may coalesce to form a bulge with 
high velocity dispersion ~ 600 Myrs into the clumpy phase 
(jElmegreen et alj 1200a ). From the kinematic sample of 13 
WiggleZ galaxies studied in Paper I, five galaxies show a sin- 
gle resolved region of emission. Two of these galaxies show 
rotation, and could be interpreted as the final coalescence 
stage predicted by the models that include clump migra- 
tion. These galaxies have average stellar masses of 3 x 10 10 
Mq and are still forming stars at a rate of 30 Mq yr~ with 
velocity dispersions of 100 km s _1 . In the kinematic sample 
3 galaxies exhibit multiple clumps, 5 are confused with ex- 
tended regions of emission and 5 show single resolved regions 
of emission. The 10 galaxies not studied in detail in this pa- 
per may also be 'clumpy' galaxies but due to projection or 
resolution effects individual regions are not distinguishable. 
However, the clumpy phase of these discs is expected to be 
short and thus we do not expect all galaxies observed to be 
in the early clumpy stage. 

A different interpretation of the high-redshift data sug- 
ests that the turbulent motions arise from star formation 



Lehnert et al.ll2009l ). This analysis, based on the correlation 
between disc velocity dispersions and Ha luminosity, is sup- 



|201Cj). We also find a tight scaling relation between lumi 
nosity and velocity dispersion (Fig. [6]) and a noisier corre- 
lation between star formation surface density and velocity 
dispersion (Fig. [SJ. However, these two parameters may be 
linked by the high dispersions within a Toomre unstable 
disc which dri ve up the star formation v ia the KS law, such 
that a oc E*' 4 IjKrumholz fc Deke]||201Cl ). as discussed in the 
previous Section. Unfortunately these arguments raise the 
problem of which is the cause and effect, high turbulence or 
high star formation rates? We note that the mechanisms we 
have suggested to generate turbulence could also be inter- 
preted as driving up star formation. 

ICeverino et alj (|201ll ) argue that because turbulent mo- 
tion can dissipate over a dynamical time, ~ 10 x shorter than 
clump lifetimes, then shortly after forming clumps must 
be supported against total free-fall collapse by internal ro- 
tation. Rotati on also provides s o me support t o clum ps in 
the models o f [Dekel et all (|2009l ). lAgertz et al.1 l|2009l ). and 
lAumer et alj (|20ld ). In simulations, clumps have rotation 
velocities ~ 120 km s _1 and v elocity dispersions co mpara- 
ble to disc velocity dispersions (|Ceverino et alj|201ll ). How- 
ever, in this simulation the rotational signature of clumps 
is found to be reduced by a factor of 7 in observations for 
beam-smearing of 0.1 arcsec, comparable to the data from 
current IFS studies. This implies that the expected observed 
velocity gradients of clumps, V ro t ~ 10 — 40 km s _1 , are not 
measurable above the errors with current instrumentation. 
The clumps extracted from discs in these simulations have 
similar properties to the observed star-forming regions, as 
shown in the first panel of Fig. [5] The simulated and ob- 
served regions overlap in size and dispersion with a similar 
scatted in velocity dispersion at a given radius. 

We test this model using the prescription of 
ICeverino et al.1 (|201ll ) which requires that clumps are rota- 
tionally supported when Kot/Kfrc > 0-5 or V T0 t/crm > \/2. 
We estimate for the WiggleZ clumps Kot/ciD in which 
uid ~ o" nct and Kot ~ «shear x 7 (correcting for beam- 
smearing assuming by a factor of 7) where Vshear is the ve- 
locity shear across the clumps as measure from the average 
of the maximum and minimum velocity residuals within the 
clump radius (~ 20 — 80 km s _1 for the WiggleZ clumps). 
Using this prescription we find that all regions could be rota- 
tionally supported. However, these conclusions are uncertain 
given our measurement errors. Observed clumps at z ~ 2 
show tentative evidence for rotational support for a fraction 
of the sample, with velocity shear measured from residual 
Ha velocity maps after model discs have been subtracted 
(jGenzel et alj|201lf ). Together, these results suggest rotation 
is not likely to be the primary support mechanism for the 
observed clumps. If clumps survive long enough to virialise 
then they could become rotationally supported. A more suit- 
able test would be to observe the kinematics of more evolved 
clumps, seen in broadband imaging, which may have had 
time to virialise. 

In conclusion, there are many processes internal and ex- 
ternal to galaxies that are able to inject energy into the ISM 
to maintain the high observed velocity dispersions. The Wig- 
gleZ clumps have consistent low metallicity across the discs 
they are embedded in, indicating that they are young with 
common evolution histories. They do not show evidence of 
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rotation suggesting that they have not virialised. However, 
this interpretation is limited by our data due to the effects 
of beam-smearing. Higher S/N and higher resolution data 
of many more clumps are needed to provide more rigourous 
tests of the existing theories of clump evolution at high red- 
shift. 

7 CONCLUSIONS 

We present 8 new clumps embedded within unstable disc 
galaxies at z ~ 1.3 from the kinematic sub-sample of the 
WiggleZ Dark Energy Survey (|Wisnioski et al.ll20liT ). The 
clumps have the following properties: 

• Average clump size of r coro ~ 1.5 kpc. 

• Average clump velocity dispersions of a ~ 90 km s , 
higher than all other clump samples at z > 1. 

• Average Jeans masses of 4.2 x 10 9 Mq, in total account- 
ing for 20-30% of the stellar mass of the discs. 

• Estimated gas mass surface densities in the range E g = 
300-2000 M Q pc~ 2 . 

• Clump star formation surface densities of E* = 0.4 — 2 
Mq yr _1 kpc~ 2 , with total clump star formation accounting 
for 30 - 40% of star formation in the discs (BG03 IMF). 

• In the stacked spectrum of all clumps we detect a broad 
Ha component of FWHM ~ 490 km s _1 which contributes 
40% to the overall clump flux and is likely due to the pres- 
ence of large-scale stellar winds driven by intense star for- 
mation. This conclusion is supported by the condition that 
all clumps are above the threshold required to drive galactic 
scale winds, E* > 0.1 M Q yr _1 kpc~ 2 . 

Clumps, HII regions, and giant HII regions are found 
to follow the same scaling relations for Ha size, velocity 
dispersion, and luminosity, indicating that they may be 
scaled-up versions of local HII regions. Our best-fitting rela- 
tions are log(cr) = (0.42 ± 0.03) x log(d) + (0.33 ± 0.09), 
log(L Ha ) = (2.72 ± 0.04) x log(d) + (31.99 ± 0.08), and 
log(L Ha ) = (4.18 ± 0.21) x log(cr) + (33.61 ± 0.31). We show 
that these relationships hold over four orders of magnitude 
in mass (10 6 — 10 10 Mq) from HII regions to giant clumps at 
high redshift. The size- velocity dispersion relation is consis- 
tent with Jeans collapse in an isothermal disc with a range in 
star formation surface density of ~ 4000 Mq pc -2 . The re- 
gions closely approximate the expected relation for idealised 
Stromgren spheres although they may not be completely 
ionisation bound as a slightly shallower relation is observed. 
This could be due to a range of factors including resolution, 
cloud density, dust, metallicity, and magnetic fields. 

The relationship between luminosity and velocity dis- 
persion is not as well understood. We find an empirical re- 
lation of L oc <7 4 , which deviates from the predicted rela- 
tion L oc a assuming that star-forming regions form un- 
der Jeans collapse and their luminosities are driven by the 
Stromgren sphere physics. A possible explanation for the 
empirical relation is that the luminosity closely traces the 
mass, such that M oc cr 4 /G 2 E, a scaling obtained for both 
the Jeans and virial masses. We find that the Jeans mass 
yields a mass in better agreement with models than the 
dynamical mass and thus is the more likely driver of the 
observed correlation. 

The data is consistent with high surface densities as a 
result of higher velocity dispersions in a marginally unstable 



disc a ssuming a Kennicutt-Schmidt law (|Krumholz fc Dekell 
|2010j). This is a reasonable approximation given that we 
find tentative evidence that star-forming regions follow the 
Kennicutt-Schmidt law. If this star-formation law holds then 
an expected range in the gas surface density would explain 
the high scatter in the a — r relation and would predict gas 
surface densities in agreement with estimates from the Jeans 
equations. 

We do not find a luminosity offset between our clumps 
and loc al HII regions as reported for other hig h-redshift 
studies (jSwinbank et all |2009| ; I Jones et al1l20ld ). We ex- 
tract measurements from the literature for giant HII re- 
gions in isolated spirals, (U)LIRGS, irregulars, and dwarf 
irregulars. After correcting historical data for updated dis- 
tances and cosmologies we find that HII regions and giant H 
II regions follow a single size-luminosity relationship. A few 
outliers persist locally from the LI RG merger of NGC4038 
and NG C4039 (iBastian et al.ll2006h and sub-L* galaxies at 
z ~ 2 Oones et al.l I2010T) and an L* galaxy at z ~ 5 
(jSwinbank et al.l I2009T ) . This offset could be due to a dif- 
ference in how regions are measured, as both high-redshift 
studies used the isophote method rather than profile fitting 
and the determination of sizes in the low-redshift sample is 
unclear. 

We determine that the most consistent method to mea- 
sure the size of an HII region is from its core radius deter- 
mined by fitting a Gaussian profile to the central brightest 
region of Ha emission. Theoretically, this measurement is 
based on the formation of idealised fully ionised regions that 
are not contaminated along the line of sight with intervening 
gas and dust. Observationally, this measurement is insensi- 
tive to variations in the background, however it becomes 
difficult to measure in low signal-to-noise data. 

The results presented here are in relative agreement 
with massive star-forming regions forming out of Toomre- 
unstable discs with high local velocity dispersion, with mea- 
sured sizes and estim ated masses in agreement with t he 
model predictions (e.g . [Noguchi 1999; Iinmc li et aill2004al lbl; 
iBournaud et al.ll2007l ; lElmegreen et alj|200a ). The data sup- 
ports models that employ radiation feedback in the form of 
large ste llar winds driven by regions of enhanced st ar for- 
mation (|Krumholz fc Dekelll2010l ; iGenel et~aI1l2010l ), how- 
ever we cannot say on what time scale the winds would 
destroy the clumps. The dominant mechanism stabilising 
clumps from rapid collapse is unclear. Given the importance 
of high velocity dispersions in the formation of large massive 
star-forming regions and minimal signature of rotation, tur- 
bulence is likely an important contributor. For the discs to 
become initially unstable the injection of energy from cold 
accretion is the probable cause. However, different mech- 
anisms may dominate throughout the lifetime of the discs, 
such as disc instabilities and fe edback from young stars, bot h 
of which are observed here (jElmegreen fc Burkertl 120101 ) . 
There are many possible mechanisms that may produce local 
velocity dispersions — inflow, winds, shocks, gravitational 
energy, star formation — further research is needed to un- 
derstand the interplay of their significance throughout the 
lifetimes of high-redshift discs. 



© 2011 RAS, MNRAS 000, [J-?? 



Scaling Relations of Star- Forming Regions 17 



ACKNOWLEDGMENTS 

We thank the referee for providing valuable comments. EW 
thanks Max Malacari and Nadine Bachmann for their valu- 
able help on measuring clump sizes and Mark Swinbank and 
Rachel Livermore for useful discussions. 

Some of the data presented herein were obtained at the 
W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, 
the University of California and the National Aeronautics 
and Space Administration. The Observatory was made pos- 
sible by the generous financial support of the W.M. Keck 
Foundation. The authors wish to recognise and acknowledge 
the very significant cultural role and reverence that the sum- 
mit of Mauna Kea has always had within the indigenous 
Hawaiian community. We are most fortunate to have the 
opportunity to conduct observations from this mountain. 



REFERENCES 

Agertz, O., Teyssier, R., & Moore, B. 2009, MNRAS, 397, 

L64 
Alexander, D. M., Swinbank, A. M., Smail, I., McDermid, 

R., & Nesvadba, N. P. H. 2010, MNRAS, 402, 2211 
Arsenault, R., & Roy, J. 1988, A&A, 201, 199 
Aumer, M., Burkert, A., Johansson, P. H., & Genzel, R. 

2010, ApJ, 719, 1230 
Baldry, I. K, & Glazebrook, K. 2003, ApJ, 593, 258 
Bastian, N., Emsellem, E., Kissler-Patig, M., & Maraston, 

C. 2006, A&A, 445, 471 
Beckman, J. E., Rozas, M., Zurita, A., Watson, R. A., & 

Knapen, J. H. 2000, AJ, 119, 2728 
Binney, J., & Tremaine, S. 2008, Galactic Dynamics: Sec- 
ond Edition, ed. Binney, J. & Tremaine, S. (Princeton 

University Press) 
Blitz, L., Fukui, Y., Kawamura, A., Leroy, A., Mizuno, N., 

& Rosolowsky, E. 2007, Protostars and Planets V, 81 
Bolatto, A. D., Leroy, A. K, Rosolowsky, E., Walter, F., 

& Blitz, L. 2008, ApJ, 686, 948 
Bournaud, F., et al. 2008, A&A, 486, 741 
Bournaud, F., Dekel, A., Teyssier, R., Cacciato, M., Daddi, 

E., Juneau, S., & Shankar, F. 2011 
Bournaud, F., Elmegreen, B. C, & Elmegreen, D. M. 2007, 

ApJ, 670, 237 
Burkert, A., et al. 2010, ApJ, 725, 2324 
Ceverino, D., Dekel, A., & Bournaud, F. 2010, MNRAS, 

404, 2151 
Ceverino, D., Dekel, A., Mandelker, N., Bournaud, F., 

Burkert, A., Genzel, R., & Primack, J. 2011, ArXiv e- 

prints 
Combes, F., Garcia-Burillo, S., Braine, J., Schinnerer, E., 

Walter, F., & Colina, L. 2011, A&A, 528, A124+ 
Cowie, L. L., Hu, E. M., & Songaila, A. 1995, AJ, 110, 1576 
Daddi, E., et al. 2010, ApJ, 713, 686 
Davies, R., et al. 2011, ApJ, 741, 69 
Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 
— . 2008, MNRAS, 383, 119 

Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785 
Dib, S., Bell, E., & Burkert, A. 2006, ApJ, 638, 797 
Dickinson, M. 2000, in Building Galaxies; from the Primor- 
dial Universe to the Present, ed. F. Hammer, T. X. Thuan, 

V. Cayatte, B. Guiderdoni, & J. T. Thanh Van , 257-+ 



Drinkwater, M. J., et al. 2010, MNRAS, 401, 1429 

Elmegreen, B. G. 1989, ApJ, 338, 178 

Elmegreen, B. G., Bournaud, F., & Elmegreen, D. M. 2008, 

ApJ, 688, 67 
Elmegreen, B. G., & Burkert, A. 2010, ApJ, 712, 294 
Elmegreen, B. G., & Elmegreen, D. M. 2005, ApJ, 627, 632 
Elmegreen, D. M., Elmegreen, B. G., Marcus, M. T., 

Shahinyan, K, Yau, A., & Petersen, M. 2009, ApJ, 701, 

306 
Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 
Ferguson, A. M. N., Wyse, R. F. C, Gallagher, III, J. S., 

& Hunter, D. A. 1996, AJ, 111, 2265 
Forster Schreiber, N. M., et al. 2009, ApJ, 706, 1364 
— . 2006, ApJ, 645, 1062 

Forster Schreiber, N. M., Shapley, A. E., Erb, D. K, Gen- 
zel, R., Steidel, C. C, Bouche, N., Cresci, G., & Davies, 

R. 2010, ArXiv e-prints 
Forster Schreiber, N. M., et al. 2011, ApJ, 739, 45 
Fuentes-Masip, O., Muhoz-Tuhon, C, Castaheda, H. O., & 

Tenorio-Tagle, G. 2000, AJ, 120, 752 
Gallagher, J. S., & Hunter, D. A. 1983, ApJ, 274, 141 
Genel, S., Bouche, N., Naab, T., Sternberg, A., & Genzel, 

R. 2010, ApJ, 719, 229 
Genel, S., et al. 2010 
Genzel, R., et al. 2011, ApJ, 733, 101 
— . 2006, Nat, 442, 786 
Green, A. "W., et al. 2010, Nat, 467, 684 
Guo, Y., Giavalisco, M., Ferguson, H. C, Cassata, P., & 

Koekemoer, A. M. 2011, ArXiv e-prints 
Gutierrez, L., Beckman, J. E., & Buenrostro, V. 2011, AJ, 

141, 113 
Heckman, T. M., Armus, L., & Miley, G. K. 1990, ApJS, 

74, 833 
Heckman, T. M., Lehnert, M. D., Strickland, D. K, & Ar- 
mus, L. 2000, ApJS, 129, 493 
Heyer, M., Krawczyk, C, Duval, J., & Jackson, J. M. 2009, 

ApJ, 699, 1092 
Immeli, A., Samland, M., Gerhard, O., & Westera, P. 

2004a, A&A, 413, 547 
Immeli, A., Samland, M., Westera, P., & Gerhard, O. 

2004b, ApJ, 611, 20 
Jones, T. A., Swinbank, A. M., Ellis, R. S., Richard, J., & 

Stark, D. P. 2010, MNRAS, 404, 1247 
Joung, M. K. R., & Mac Low, M.-M. 2006, ApJ, 653, 1266 
Joung, M. R., Mac Low, M.-M., & Bryan, G. L. 2009, ApJ, 

704, 137 
Kennicutt, Jr., R. C. 1979, ApJ, 228, 696 
— . 1998, ARA&A, 36, 189 

Kennicutt, Jr., R. C, et al. 2003, PASP, 115, 928 
— . 2007, ApJ, 671, 333 
Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, 

MNRAS, 363, 2 
Krumholz, M., & Burkert, A. 2010, ApJ, 724, 895 
Krumholz, M. R., & Dekel, A. 2010, MNRAS, 406, 112 
Larkin, J., et al. 2006, New A Rev., 50, 362 
Larson, R. B. 1981, MNRAS, 194, 809 
Law, D. R., Steidel, C. C, Erb, D. K, Larkin, J. E., Pettini, 

M., Shapley, A. E., & Wright, S. A. 2007, ApJ, 669, 929 
Le Than, L., Lehnert, M. D., Di Matteo, P., Nesvadba, 

N. P. H., & van Driel, W. 2011, A&A, 530, L6 
Lehnert, M. D., Nesvadba, N. P. H., Than, L. L., Di Mat- 



© 2011 RAS, MNRAS 000, [TV?? 



18 E. Wisnioski et al. 

teo, P., van Driel, W., Douglas, L. S., Chemin, L., & Bour- 

naud, F. 2009, ApJ, 699, 1660 
Mac Low, M.-M., & Klessen, R. S. 2004, Reviews of Modern 

Physics, 76, 125 
Mac Low, M.-M., Klessen, R. S., Burkert, A., & Smith, 

M. D. 1998, Physical Review Letters, 80, 2754 
Melnick, J., Moles, M., Terlevich, R., & Garcia-Pelayo, J.- 

M. 1987, MNRAS, 226, 849 
Monreal-Ibero, A., Colina, L., Arribas, S., & Garcia-Marm, 

M. 2007, A&A, 472, 421 
Murray, N., Quataert, E., & Thompson, T. A. 2010, ApJ, 

709, 191 
Noguchi, M. 1999, ApJ, 514, 77 
Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae 

and active galactic nuclei, ed. Osterbrock, D. E. 
Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Pleuss, P. O., Heller, C. H., & Fricke, K. J. 2000, A&A, 

361, 913 
Relano, M., Beckman, J. E., Zurita, A., Rozas, M., & Gi- 

ammanco, C. 2005, A&A, 431, 235 
Rozas, M., Richer, M. G., Lopez, J. A., Relano, M., & 

Beckman, J. E. 2006, A&A, 455, 539 
Samland, M., & Gerhard, O. E. 2003, A&A, 399, 961 
Sandage, A., & Tammann, G. A. 1974, ApJ, 190, 525 
Schmidt, M. 1959, ApJ, 129, 243 
Shapiro, K. L., et al. 2009, ApJ, 701, 955 
Smith, M. G., & Weedman, D. W. 1970, ApJ, 161, 33 
Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 

1987, ApJ, 319, 730 
Spitzer, Jr., L. 1942, ApJ, 95, 329 
Stone, J. M., Ostriker, E. C, & Gammie, C. F. 1998, ApJL, 

508, L99 
Swinbank, A. M., et al. 2010, Nat, 464, 733 
— . 2009, MNRAS, 400, 1121 
Tacconi, L. J., et al. 2010, Nat, 463, 781 
Tenorio-Tagle, G., Munoz-Tunon, C, & Cox, D. P. 1993, 

ApJ, 418, 767 
Terlevich, R., & Melnick, J. 1981, MNRAS, 195, 839 
Tiran, L. L., Lehnert, M. D., van Driel, W., Nesvadba, N. 

P. H., & Matteo, P. D. 2011 
Toomre, A. 1964, ApJ, 139, 1217 
Tully, R. B., Rizzi, L., Shaya, E. J., Courtois, H. M., 

Makarov, D. I., & Jacobs, B. A. 2009, AJ, 138, 323 
van Dam, M. A., et al. 2006, PASP, 118, 310 
van den Bergh, S., Abraham, R. G., Ellis, R. S., Tanvir, 

N. R., Santiago, B. X., & Glazebrook, K. G. 1996, AJ, 

112, 359 
van Starkenburg, L., van der Werf, P. P., Franx, M., Labbe, 

I., Rudnick, G., & Wuyts, S. 2008, A&A, 488, 99 
Wada, K., Meurer, G., & Norman, C. A. 2002, ApJ, 577, 

197 
Wisnioski, E., et al. 2011, MNRAS, 417, 2601 
Wizinowich, P. L., et al. 2006, PASP, 118, 297 



© 2011 RAS, MNRAS 000, [TL?? 



Scaling Relations of Star-Forming Regions 19 



APPENDIX A: CORRECTIONS TO CLUMP & 
HII REGION MEASUREMENTS 

This appendix summarises the methods employed to 
calculate sizes, luminosities, and velocity dispersions from 
the literature, used for comparison with WiggleZ clumps. 
Details of the corrections we introduced to improve the 
accuracy of comparison are included. 



Al HII Regions 

SINGS measurements 

From the SINGS survey we measure HII region size and 
luminosities from NGC 24, NGC 628, NGC 925, NGC 1566, 
NGC 4254, NGC 3938, NGC 7552 as detailed in Section 3. 
These galaxies were chosen for the availability of high quality 
Ha images and high numbers of HII regions. 



A2 Giant HII regions 

Gallagher & Hunter 1983 

The Ha flux is measured within a fixed 200 pc radius for 
each HII region. This essentially yields the flux within 
the core radius. The authors estimate that there is up to 
a factor of two error in the flux measurements. Velocity 
dispersions are corrected for Doppler broadening and 
instrumental broadening. We recalculate the luminosity 
and diameters with Ho = 70 km s" 1 Mpc" 1 from Ho — 50 
km s" 1 Mpc -1 originally used in the paper, using new 
distances from the Extragalactic Distance Database. 

Arsenault 1988 

Recalculation of HII region luminosities and diameters 
with Ho = 70 km s~ Mpc - and new distances from the 
Extragalactic Distance Database is readily possible as the 
authors provide the original measurements of angular di- 
ameter and Ha flux. Velocity dispersions are not measured 
for this sample. 

Bastian et al. 2006 

Six HII regions are observed in the Antennae galaxies, NGC 
4038 and NGC 4039, with the VLT-VIMOS IFS. The H 7 , 
H/3, and [OIII] emission lines are observed. The Antennae 
is classified as a LIRG. Luminosities are inferred from 
H/3 and converted to Ha assuming case B recombination 
(10 000 K). We recalculate the luminosity and diameters 
with Ho = 70 km s" 1 Mpc" 1 from H = 75 km s" 1 
Mpc - originally used in the paper, with a new distance 
of D = 19.2 Mpc from the Extragalactic Distance Database. 

Rozas et al. 2006 

This sample includes only the intrinsically brightest H 
II regions within NGC 157, NGC 925, NGC 3631, NGC 
6764, NGC 3344, NGC 4321, NGC 5364, NGC 5055, NGC 
5985, and NGC 7479. The purpose of the selection was to 
choose the youngest, most massive normal star clusters to 
minimise the processes available to drive the HII region 
away from virial equilibrium. To find an accurate velocity 
dispersion the authors determine the optimal number of 
Gaussian components to fit to each line profile using the 
component which contributes the most luminosity in each H 



II region. For region size they adopt one half t he total radius 
iven in HII region catalogues as specified bv lRelano et al.l 
2005! ) , which is based on the 40% isophote method of 
McCall 1990, which is equivalent to our core method. The 
velocity dispersion is measured within this radius and was 
corrected for natural broadening, Doppler broadening, and 
instrumental broadening. We recalculate the luminosity and 
diameters with Ho = 70 km s" 1 Mpc -1 from Ho = 75 km 
s" 1 Mpc -1 originally used in the paper, with new distances 
from the Extragalactic Distance Database. 

Monreal-Ibero et al. 2007 

The regions selected from this paper are tidal dwarf galaxy 
candidates. However, they are selected in a similar manner 
to other high-redshift and low-redshift regions used for com- 
parison in this paper and as such we include them as giant 
HII regions. They are selected within 5 local ULIRGs with 
WFPC2 and IFS data. Isophotal and effective radii are mea- 
sured for all regions. We use the effective radii as it better 
approximates the core method. Luminosity is recalculated 
from the given flux values without extinction correction. 



A3 High-Redshift Clumps 

SINS Survey: Genzel et al. (2011) &z Forster- 
Schrieber et al. (2011) 

The clump data taken from iGenzel et al.l (|201ll ) are for 
BX482-clumpA, ZC782941-clumpA, ZC406690-clumpA, 
ZC406690-clumpB, and ZC406690-clumpC. Integrated 
clump properties such as BX599-all and D3al5504- 
clumpsA-F are excluded. Clump radii are given by the 
HWHM from Gaussian fits to clumps in velocity channels 
of the IFS data cube. The instrumental resolution is sub- 
tracted in quadrature. We convert the HWHM to the r.m.s. 
of a Gaussian profile, the radius used for the core method, 
by r CO rc = 2 ?y H y . The velocity dispersion of the clumps 
is the intrinsic local velocity dispersion after removal of 

beam-smeared rotation and instru mental resolution. 

The clump data taken from IForster Schreiber et al.l 

(|2011f ) are exclusively from BX482. Clump-1 in this sample 
is defined to be th e sam e clump as BX482-clu m pA in 
I Genzel et all (|201ll ). In IForster Schreiber et all (|201ll ) 
the flux, F ohs = 3.92 ± 0.05 x 10" 17 ergs s" 1 cm" 2 is 
measur ed within r P ho t = 2. 3 kpc and has -Rhwhm = 0.96 
kpc in I Genzel et al.l (|2011f ) the -Rhwhm = 1 kpc and 
F hs = 3.5 x 10" 17 ergs s" 1 cm" 2 . For comparison purposes 



we use r ph ot 



radii. 



Swinbank et a l. 2009 &z Jo nes et al. 2 010 

ISwinbank et all <|2009h and Ijones et all (|2010T ) publish 
clumps in L* lensed galaxies at z ~ 2 — 4 from IFS data. 
They derive sizes using the isophote method. We correct 
the cosmology to Ho — 70 km s" 1 Mpc" 1 and Qm = 0.27, 
Ov = 0.73. 
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